Balancing economic and ecological functions in smallholder and industrial oil palm plantations.

The expansion of the oil palm industry in Indonesia has improved livelihoods in rural communities, but comes at the cost of biodiversity and ecosystem degradation. Here, we investigated ways to balance ecological and economic outcomes of oil palm cultivation. We compared a wide range of production systems, including smallholder plantations, industrialized company estates, estates with improved agronomic management, and estates with native tree enrichment. Across all management types, we assessed multiple indicators of biodiversity, ecosystem functions, management, and landscape structure to identify factors that facilitate economic-ecological win-wins, using palm yields as measure of economic performance. Although, we found that yields in industrialized estates were, on average, twice as high as those in smallholder plantations, ecological indicators displayed substantial variability across systems, regardless of yield variations, highlighting potential for economic-ecological win-wins. Reducing management intensity (e.g., mechanical weeding instead of herbicide application) did not lower yields but improved ecological outcomes at moderate costs, making it a potential measure for balancing economic and ecological demands. Additionally, maintaining forest cover in the landscape generally enhanced local biodiversity and ecosystem functioning within plantations. Enriching plantations with native trees is also a promising strategy to increase ecological value without reducing productivity. Overall, we recommend closing yield gaps in smallholder cultivation through careful intensification, whereas conventional plantations could reduce management intensity without sacrificing yield. Our study highlights various pathways to reconcile the economics and ecology of palm oil production and identifies management practices for a more sustainable future of oil palm cultivation.

Decoding autumn phenology: Unraveling the link between observation methods and detected environmental cues.

Leaf coloring and fall mark the end of the growing season (EOS), playing essential roles in nutrient cycling, resource allocation, ecological interactions, and as climate change indicators. However, understanding future changes in autumn phenology is challenging due to the multitude of likely environmental cues and substantial variations in timing caused by different derivation methods. Yet, it remains unclear whether these two factors are independent or if methodological uncertainties influence the environmental cues determined. We derived start of growing season (SOS) and EOS at a mixed beech forest in Central Germany for the period 2000-2020 based on four different derivation methods using a unique long-term data set of in-situ data, canopy imagery, eddy covariance measurements, and satellite remote sensing data and determined their influence on a predictor analysis of leaf senescence. Both SOS and EOS exhibited substantial ranges in mean onset dates (39.5 and 28.6 days, respectively) across the different methods, although inter-annual variations and advancing SOS trends were similar across methods. Depending on the data, EOS trends were advanced or delayed, but inter-annual patterns correlated well (mean r = .46). Overall, warm, dry, and less photosynthetically productive growing seasons were more likely to be associated with a delayed EOS, while colder, wetter, and more photosynthetically productive vegetation periods resulted in an earlier EOS. In addition, contrary to recent results, no clear influence of pre-solstice vegetation activity on the timing of senescence was detected. However, most notable were the large differences in sign and strength of potential drivers both in the univariate and multivariate analyses when comparing derivation methodologies. The results suggest that an ensemble analysis of all available phenological data sources and derivation methods is needed for general statements on autumn phenology and its influencing variables and correct implementation of the senescence process in ecosystem models.

Transformation scenarios towards multifunctional landscapes: A multi-criteria land-use allocation model applied to Jambi Province, Indonesia.

In tropical regions, shifting from forests and traditional agroforestry to intensive plantations generates conflicts between human welfare (farmers' demands and societal needs) and environmental protection. Achieving sustainability in this transformation will inevitably involve trade-offs between multiple ecological and socioeconomic functions. To address these trade-offs, our study used a new methodological approach allowing the identification of transformation scenarios, including theoretical landscape compositions that satisfy multiple ecological functions (i.e., structural complexity, microclimatic conditions, organic carbon in plant biomass, soil organic carbon and nutrient leaching losses), and farmers needs (i.e., labor and input requirements, total income to land, and return to land and labor) while accounting for the uncertain provision of these functions and having an actual potential for adoption by farmers. We combined a robust, multi-objective optimization approach with an iterative search algorithm allowing the identification of ecological and socioeconomic functions that best explain current land-use decisions. The model then optimized the theoretical land-use composition that satisfied multiple ecological and socioeconomic functions. Between these ends, we simulated transformation scenarios reflecting the transition from current land-use composition towards a normative multifunctional optimum. These transformation scenarios involve increasing the number of optimized socioeconomic or ecological functions, leading to higher functional richness (i.e., number of functions). We applied this method to smallholder farms in the Jambi Province, Indonesia, where traditional rubber agroforestry, rubber plantations, and oil palm plantations are the main land-use systems. Given the currently practiced land-use systems, our study revealed short-term returns to land as the principal factor in explaining current land-use decisions. Fostering an alternative composition that satisfies additional socioeconomic functions would require minor changes ("low-hanging fruits"). However, satisfying even a single ecological indicator (e.g., reduction of nutrient leaching losses) would demand substantial changes in the current land-use composition ("moonshot"). This would inevitably lead to a profit decline, underscoring the need for incentives if the societal goal is to establish multifunctional agricultural landscapes. With many oil palm plantations nearing the end of their production cycles in the Jambi province, there is a unique window of opportunity to transform agricultural landscapes.

Drought alters aboveground biomass production efficiency: Insights from two European beech forests.

The fraction of photosynthetically assimilated carbon that trees allocate to long-lasting woody biomass pools (biomass production efficiency - BPE), is a key metric of the forest carbon balance. Its apparent simplicity belies the complex interplay between underlying processes of photosynthesis, respiration, litter and fruit production, and tree growth that respond differently to climate variability. Whereas the magnitude of BPE has been routinely quantified in ecological studies, its temporal dynamics and responses to extreme events such as drought remain less well understood. Here, we combine long-term records of aboveground carbon increment (ACI) obtained from tree rings with stand-level gross primary productivity (GPP) from eddy covariance (EC) records to empirically quantify aboveground BPE (= ACI/GPP) and its interannual variability in two European beech forests (Hainich, DE-Hai, Germany; Sorø, DK-Sor, Denmark). We found significant negative correlations between BPE and a daily-resolved drought index at both sites, indicating that woody growth is de-prioritized under water limitation. During identified extreme years, early-season drought reduced same-year BPE by 29 % (Hainich, 2011), 31 % (Sorø, 2006), and 14 % (Sorø, 2013). By contrast, the 2003 late-summer drought resulted in a 17 % reduction of post-drought year BPE at Hainich. Across the entire EC period, the daily-to-seasonal drought response of BPE resembled that of ACI, rather than that of GPP. This indicates that BPE follows sink dynamics more closely than source dynamics, which appear to be decoupled given the distinctive climate response patterns of GPP and ACI. Based on our observations, we caution against estimating the magnitude and variability of the carbon sink in European beech (and likely other temperate forests) based on carbon fluxes alone. We also encourage comparable studies at other long-term EC measurement sites from different ecosystems to further constrain the BPE response to rare climatic events.

Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown.

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Leaf-level coordination principles propagate to the ecosystem scale.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories - the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis - are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

Differences between tree stem CO2 efflux and O2 influx rates cannot be explained by internal CO2 transport or storage in large beech trees.

Tree stem respiration (RS ) is a substantial component of the forest carbon balance. The mass balance approach uses stem CO2 efflux and internal xylem fluxes to sum up RS , while the oxygen-based method assumes O2 influx as a proxy of RS . So far, both approaches have yielded inconsistent results regarding the fate of respired CO2 in tree stems, a major challenge for quantifying forest carbon dynamics. We collected a data set of CO2 efflux, O2 influx, xylem CO2 concentration, sap flow, sap pH, stem temperature, nonstructural carbohydrates concentration and potential phosphoenolpyruvate carboxylase (PEPC) capacity on mature beech trees to identify the sources of differences between approaches. The ratio of CO2 efflux to O2 influx was consistently below unity (0.7) along a 3-m vertical gradient, but internal fluxes did not bridge the gap between influx and efflux, nor did we find evidence for changes in respiratory substrate use. PEPC capacity was comparable with that previously reported in green current-year twigs. Although we could not reconcile differences between approaches, results shed light on the uncertain fate of CO2 respired by parenchyma cells across the sapwood. Unexpected high values of PEPC capacity highlight its potential relevance as a mechanism of local CO2 removal, which merits further research.

Joint optimization of land carbon uptake and albedo can help achieve moderate instantaneous and long-term cooling effects.

Both carbon dioxide uptake and albedo of the land surface affect global climate. However, climate change mitigation by increasing carbon uptake can cause a warming trade-off by decreasing albedo, with most research focusing on afforestation and its interaction with snow. Here, we present carbon uptake and albedo observations from 176 globally distributed flux stations. We demonstrate a gradual decline in maximum achievable annual albedo as carbon uptake increases, even within subgroups of non-forest and snow-free ecosystems. Based on a paired-site permutation approach, we quantify the likely impact of land use on carbon uptake and albedo. Shifting to the maximum attainable carbon uptake at each site would likely cause moderate net global warming for the first approximately 20 years, followed by a strong cooling effect. A balanced policy co-optimizing carbon uptake and albedo is possible that avoids warming on any timescale, but results in a weaker long-term cooling effect.

Land-use change shifts and magnifies seasonal variations of the decomposer system in lowland tropical landscapes.

Deforestation and agricultural expansion in the tropics affect local and regional climatic conditions, leading to synergistic negative impacts on land ecosystems. Climatic changes manifest in increased inter- and intraseasonal variations and frequency of extreme climatic events (i.e., droughts and floods), which have evident consequences for aboveground biodiversity. However, until today, there have been no studies on how land use affects seasonal variations below ground in tropical ecosystems, which may be more buffered against climatic variation. Here, we analyzed seasonal variations in soil parameters, basal respiration, microbial communities, and abundances of soil invertebrates along with microclimatic conditions in rainforest and monocultures of oil palm and rubber in Sumatra, Indonesia. About 75% (20 out of 26) of the measured litter and soil, microbial, and animal parameters varied with season, with seasonal changes in 50% of the parameters depending on land use. Land use affected seasonal variations in microbial indicators associated with carbon availability and cycling rate. The magnitude of seasonal variations in microbial parameters in the soil of monocultures was almost 40% higher than in the soil of rainforest. Measured parameters were associated with short-term climatic conditions (3-day period air humidity) in plantations, but not in rainforest, confirming a reduced soil buffering ability in plantations. Overall, our findings suggest that land use temporally shifts and increases the magnitude of seasonal variations of the belowground ecosystem compartment, with microbial communities responding most strongly. The increased seasonal variations in soil biota in plantations likely translate into more pronounced fluctuations in essential ecosystem functions such as nutrient cycling and carbon sequestration, and these ramifications ultimately may compromise the stability of tropical ecosystems in the long term. As the observed seasonal dynamics is likely to increase with both local and global climate change, these shifts need closer attention for the long-term sustainable management of plantation systems in the tropics.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems.

Understanding the critical soil moisture (SM) threshold (θcrit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θcrit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θcrit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θcrit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θcrit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θcrit . We further found that soil texture dominates the spatial variability of θcrit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θcrit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

Dataset on microclimate and drone-based thermal patterns within an oil palm agroforestry system.

Microclimate and Land Surface Temperature (LST) are important analytical variables used to understand complex oil palm agroforestry systems and their effects on biodiversity and ecosystem functions. In order to examine experimental effects of tree species richness (0, 1, 2, 3 or 6), plot size (25 m2, 100 m2, 400 m2, 1600 m2) and stand structural complexity on microclimate and Land Surface Temperature, related data were collected following a strict design. The experiment was carried out in the Jambi province, in Sumatra (Indonesia), as part of the collaborative project EFForTS [Ecological and Socioeconomic Functions of Tropical Lowland Rainforest Transformation Systems]. Microclimate data collected using miniaturized data loggers combined with drone-based thermal data were considered within an oil palm plantation enriched with six target tree species. The timeframe considered for data analysis was 20th September 2017 to 26th September 2017. The experiment data can be used for comparison with data from conventional oil palm agroforestry systems in the tropics. They can more specifically be used as reference to assess microclimate and Land Surface Temperature patterns within similar agroforestry systems.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

Sensitivity of gross primary productivity to climatic drivers during the summer drought of 2018 in Europe.

In summer 2018, Europe experienced a record drought, but it remains unknown how the drought affected ecosystem carbon dynamics. Using observations from 34 eddy covariance sites in different biomes across Europe, we studied the sensitivity of gross primary productivity (GPP) to environmental drivers during the summer drought of 2018 versus the reference summer of 2016. We found a greater drought-induced decline of summer GPP in grasslands (-38%) than in forests (-10%), which coincided with reduced evapotranspiration and soil water content (SWC). As compared to the 'normal year' of 2016, GPP in different ecosystems exhibited more negative sensitivity to summer air temperature (Ta) but stronger positive sensitivity to SWC during summer drought in 2018, that is, a stronger reduction of GPP with soil moisture deficit. We found larger negative effects of Ta and vapour pressure deficit (VPD) but a lower positive effect of photosynthetic photon flux density on GPP in 2018 compared to 2016, which contributed to reduced summer GPP in 2018. Our results demonstrate that high temperature-induced increases in VPD and decreases in SWC aggravated drought impacts on GPP. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018.

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO2 exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO2 uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Trade-offs between multifunctionality and profit in tropical smallholder landscapes.

Land-use transitions can enhance the livelihoods of smallholder farmers but potential economic-ecological trade-offs remain poorly understood. Here, we present an interdisciplinary study of the environmental, social and economic consequences of land-use transitions in a tropical smallholder landscape on Sumatra, Indonesia. We find widespread biodiversity-profit trade-offs resulting from land-use transitions from forest and agroforestry systems to rubber and oil palm monocultures, for 26,894 aboveground and belowground species and whole-ecosystem multidiversity. Despite variation between ecosystem functions, profit gains come at the expense of ecosystem multifunctionality, indicating far-reaching ecosystem deterioration. We identify landscape compositions that can mitigate trade-offs under optimal land-use allocation but also show that intensive monocultures always lead to higher profits. These findings suggest that, to reduce losses in biodiversity and ecosystem functioning, changes in economic incentive structures through well-designed policies are urgently needed.

Measured greenhouse gas budgets challenge emission savings from palm-oil biodiesel.

The potential of palm-oil biofuels to reduce greenhouse gas (GHG) emissions compared with fossil fuels is increasingly questioned. So far, no measurement-based GHG budgets were available, and plantation age was ignored in Life Cycle Analyses (LCA). Here, we conduct LCA based on measured CO2, CH4 and N2O fluxes in young and mature Indonesian oil palm plantations. CO2 dominates the on-site GHG budgets. The young plantation is a carbon source (1012 ± 51 gC m-2 yr-1), the mature plantation a sink (-754 ± 38 gC m-2 yr-1). LCA considering the measured fluxes shows higher GHG emissions for palm-oil biodiesel than traditional LCA assuming carbon neutrality. Plantation rotation-cycle extension and earlier-yielding varieties potentially decrease GHG emissions. Due to the high emissions associated with forest conversion to oil palm, our results indicate that only biodiesel from second rotation-cycle plantations or plantations established on degraded land has the potential for pronounced GHG emission savings.

It is not just a 'trade-off': indications for sink- and source-limitation to vegetative and regenerative growth in an old-growth beech forest.

Controls on tree growth are key issues in plant physiology. The hypothesis of our study was that the interannual variability of wood and fruit production are primarily controlled directly by weather conditions (sink limitation), while carbon assimilation (source limitation) plays a secondary role. We analyzed the interannual variability of weather conditions, gross primary productivity (GPP) and net primary productivity (NPP) of wood and fruits of an old-growth, unmanaged Fagus sylvatica forest over 14 yr, including six mast years. In a multiple linear regression model, c. 71% of the annual variation in wood-NPP could be explained by mean air temperature in May, precipitation from April to May (positive influence) and fruit-NPP (negative influence). GPP of June to July solely explained c. 42% of the variation in wood-NPP. Fruit-NPP was positively related to summer precipitation 2 yr before (R2  = 0.85), and negatively to precipitation in May (R2  = 0.83) in the fruit years. GPP had no influence on fruit-NPP. Our results suggest a complex system of sink and source limitations to tree growth driven by weather conditions and going beyond a simple carbon-mediated 'trade-off' between regenerative and vegetative growth.

Memory effects of climate and vegetation affecting net ecosystem CO2 fluxes in global forests.

Forests play a crucial role in the global carbon (C) cycle by storing and sequestering a substantial amount of C in the terrestrial biosphere. Due to temporal dynamics in climate and vegetation activity, there are significant regional variations in carbon dioxide (CO2) fluxes between the biosphere and atmosphere in forests that are affecting the global C cycle. Current forest CO2 flux dynamics are controlled by instantaneous climate, soil, and vegetation conditions, which carry legacy effects from disturbances and extreme climate events. Our level of understanding from the legacies of these processes on net CO2 fluxes is still limited due to their complexities and their long-term effects. Here, we combined remote sensing, climate, and eddy-covariance flux data to study net ecosystem CO2 exchange (NEE) at 185 forest sites globally. Instead of commonly used non-dynamic statistical methods, we employed a type of recurrent neural network (RNN), called Long Short-Term Memory network (LSTM) that captures information from the vegetation and climate's temporal dynamics. The resulting data-driven model integrates interannual and seasonal variations of climate and vegetation by using Landsat and climate data at each site. The presented LSTM algorithm was able to effectively describe the overall seasonal variability (Nash-Sutcliffe efficiency, NSE = 0.66) and across-site (NSE = 0.42) variations in NEE, while it had less success in predicting specific seasonal and interannual anomalies (NSE = 0.07). This analysis demonstrated that an LSTM approach with embedded climate and vegetation memory effects outperformed a non-dynamic statistical model (i.e. Random Forest) for estimating NEE. Additionally, it is shown that the vegetation mean seasonal cycle embeds most of the information content to realistically explain the spatial and seasonal variations in NEE. These findings show the relevance of capturing memory effects from both climate and vegetation in quantifying spatio-temporal variations in forest NEE.